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We have treated the adiabatic matter creation process in various three-brane models by 
applying thermodynamics of open systems. The matter creation rate is found to affect 
the evolution of scale factor and energy density of the universe. We find modification at 
early stages of cosmic dynamics. In GB and RS brane worlds, by chosing appropriate 
parameters we obtain standard scenario, while the warped DGP model has different 
solutions. During later stages, since the matter creation is negligible the evolution reduces 
to FRW expansion, in RS and GB models. 
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1. Matter Creation Scenario 

Einstein's field equations are adiabatic and reversible and hence cannot provide an 
explanation for the origin of cosmological entropy. 1 But matter constituents may 
be produced quantum mechanically in the framework of Einstein's equations. 2,3 
The quantum fluctuations of the gravitational field near singularity, as well as the 
dynamics of any quantum field living on such a non stationary background, lead to 
matter- creating phenomena, which can in turn modify the behavior of the early 
universe. 4 A detailed microscopic description of matter creation phenomena should 
be provided by a fully self-consistent quantum field theory on a curved space time. 5 
Many attempts have been made to treat the matter creation process at a phe- 
nomenological macroscopic level. 6-13 The basic idea is that matter creation con- 
tributes at the level of the Einstein field equations as a negative pressure term. 
Prigogine et al. incorporated its effect in the classical Einstein field equations. 1 At 
the expense of the gravitational field, matter creation can occur only as an irre- 
versible process. The additional pressure term due to the creation depends on the 
creation rate. The total entropy within a comoving volume is not a constant, due 
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to the variation of number of particles. This approach explained the average effect 
induced on the dynamics of an isotropic expanding universe. In this work we study 
the effect of adiabatic matter creation on brane world models. We haven't consid- 
ered the details of the source of created matter. Recently Tetradis has shown that 
absorption of energy flux from the bulk by the three brane leads to acceleration or 
deceleration of the universe. 14 We have also derived a phantom field equivalent for- 
mulation. This paper is organised as follows. In section 2 we give a brief description 
of RSII brane world model. In the section 3, thermodynamics of matter creation is 
explained. In section 4 we describe the effect of matter creation on brane worlds in 
general, and on RSII in section 5. An equivalent phantom field formulation is given 
in section 6. In sections 7 and 8 we give the adiabatic matter creation effects on 
DGP and GB brane models. 



2. RSII Brane Model 

The idea that our 4-dimensional universe might be a 3-brane embedded in a higher 
dimensional space-time has attracted much attention over the last several years. 
Randall and Sundrum (RS) models 15 ' 16 are the most important. According to the 
brane world scenario, the physical fields in our 4-dimensional space-time are confined 
to the three brane, while gravity can freely propagate in the bulk space-time. We 
assume a Friedmann metric for 4-dimensional spacetime embedded in the bulk 5- 
dimcnsional space with a cosmological constant A5. 
The Friedmann equation on the RSII brane is, 17,18 

with 

2 8ttG A 4 <j 2 A 5 

K d = „ , ana — = — + — , (2) 

4 3 3 36 6 ' 

where A4 is the cosmological constant on the brane, a is the brane tension, p is 
energy density of ordinary cosmological matter on the brane, k is the curvature 
index, and the term with C (a constant) is the dark radiation term. 19 Cosmology 
constrains the amount of dark radiation to be at most 10 percent of the energy 
density in photons. 20 In the following we will assume AdS bulk spacetime, and 
C = 0, k = 0, and A 4 = 

In RSII model the important change in Friedmann's equation compared to the 
usual four dimensional form is the appearance of a term proportional to p. It tells 
us that if the matter energy density is much larger than the brane tension, i.e. 
p >> a, the expansion rate H is proportional p. Only in the limit where the brane 
tension is much larger than the matter energy density, the usual behavior H oc ^fp 
is recovered. 
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3. Thermodynamics of Matter Creation in Cosmology 

The Einstein's field equations with Bianchi identities give, for homogeneous and 
isotropic universes in closed systems, 

d{pV) = -p dV. (3) 

where ft is the true thermodynamical pressure. 

This equation is used to describe an adiabatic evolution of a closed system with 
a comoving volume V . But in the presence of matter creation, the analysis is in 
the context of open systems. 1,21 In this case the number of particles N in a given 
volume is not fixed. For adiabatic transformation, 

d{pV) +pdV - ^ d{nV) = 0. (4) 

where n = y and h = p + p is the enthalpy per unit volume. 

In such a transformation, the heat received by the system is due to the change 
of the number of particles. In our cosmological context, this change is due to the 
transfer of energy from gravitation to matter. Correspondingly the entropy change 
dS, for adiabatic transformation in open systems is, 

TdS=- dinV) - n ( P V) = T— d(nV) (5) 
n n 

where p = h — T s is the chemical potential and s = y. Therefore the only particle 
number variations admitted are such that 

dN = d(nV) > (6) 

The Eq. (5) can be expressed as, 

p=-n (7) 

p= HA^ (8) 
n 

where the over dot represents time derivative. The creation of matter corresponds 
to a supplementary pressure p c , which is negative or zero depending on the presence 
or absence of particle production. Thus Eq. (5) can be written as, 

d(pV) = -( P + p c )dV = -ftdV (9) 

where p is the true thermodynamical pressure and 

h djnV) p + p d(nV) 
Pc ~ n dV ~ n dV [W) 
This scenario is applied to cosmology as follows. In the case of an isotropic and 
homogeneous universe, we choose, for V, the value V — R 3 (r) where R is the scale 
factor. The above expression becomes, 

Pc = -^£(n + 3nH) (11) 
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The conservation equation becomes, 

p = -3H(p + p + Pc ) = -3H(p + p) + 3H (n + 3nH) (12) 

11 

P=-(P + P) (13) 



with 

S N h + 3nH 



(14) 



S N n 
The balance equation for particle number density is, 1 

N l =nu l , N? a =n + 3nH = riT (15) 

with 

N = na 3 , and T = — = ^ (16) 

where p is thermostatic pressure, u l is fluid 4-velocity, n is particle number density, 
and T is matter creation rate. For T > 0, there is particle production, and T < 
corresponds to annihilation. 

The above-mentioned macroscopic phenomenological approach is not completely 
self consistent since it cannot determine the expression for the rate of particle cre- 
ation, which is an open degree of the scheme. 



4. Effect of Mater creation on Brane Worlds 

The existence of some kind of negative-pressure dark energy is usually assumed to 
explain the cosmic acceleration. Vacuum energy or a cosmological constant is the 
simplest dark energy candidate. It has been known that matter creation may lead 
to negative pressures. 22 The modified Friedmann equation on RSII brane world 
indicates that the universe evolves slowly in the very early times. The process of 
matter creation may add to the accelerated expansion and thus can compensate to 
the effect of brane world at early times. A scenario with particle production driving 
cosmic acceleration was proposed by Zimdahl et al. 23 Using a phenomenological 
macroscopic description of matter creation, 1 we investigate this scenario in the brane 
world context. However, we do not consider that particle creation as the only source 
for cosmic acceleration. It has been shown based on the supernova data that matter 
creation alone does not drive cosmic acceleration. 24 Our motivation is the following: 
the modified Friedmann equation on RSII brane world shows that there is transition 
for the universe evolution between an early high energy regime characterized by 
the behavior R oc fV 3 ( 1+w ) and a late low energy regime characterized by the 
standard evolution R cx t 2 / 3 ( 1 +") 25,18 ^ e matter crea ti n process due to its 
negative pressure can be the cause of this transition. Our attempt is to consider the 
possible effects on dynamics of the early universe due the matter creation. 
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5. Effect of Adiabatic Matter Creation on the RSII Brane Model 

The field equations on the RS II brane model are, 26 

TT o R 2 k 2 p A C k 



R 2 3 2cr 3 R 4 R 2 . 

R R 2 



R R 2 k 2 , . / p\ 2C A 



where k 2 = 87rG , and lo = 7 — 1 = 2 for perfect fluid. With particle creation, 
creation pressure, Eq. (11), is added to thermostatic pressure p. 

For simplicity we take Lambda = 0. Also, since the dark radiation term decays 
rapidly we can put C = 0, Eq. (18) becomes, 



i? i? 2 K 2 ^ ; p + }3 C 

P 

The energy conservation equation is unchanged in 5-d setup, which including 
particle creation, becomes, 



p + 'iH (p + p + Pc) = (21) 
For perfect fluid, using Eq. (19) we get, 

P + ^p(l-^] =0 (22) 



We take the linear dependence of matter creation rate on Hubble parameter H 
as, 27 

iP = (33nH (23) 

where the parameter (3 is defined in the interval [0, 1] and assumed to be constant. 
The observational value of H to (0.86 < H t n < 1.91) can be translated to the value 
of 13. The joint result of SNe la and GRB indicate that (3 is (0.502 ± 0.038). The 
GRB data shows (3 to be greater than 1/3 (ie.,0.349± 0.185), and the SNIa data 
alone shows (0.537 ± 0.040). 28 

This linear relation implies -0 cx p, at high energy densities. (Recently Freaza et. 
al. 2 have considered a general form for the matter creation in non-brane scenario 
as, ^ = 3(3nH (H/H ) a , where a and (3 are dimensionless constants). Now the 
solution of Eq. (22) is, 

p = po R~ q where q = 37 (1 - (3) (24) 

Substituting for p, and for C = k = A = 0, Eq. (17) becomes, 
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Denning R q — X, this becomes, 

A „ „1 ^2 Jl „i 



X 2 = K Poq X+ K Poq =AX + B (26) 
3 6 <7 



t-t = ^X + l (27) 



Solving we get 



Using the boundary condition, X = R q = 0, when t — , we find the solution for 
the scale factor as follows, 

R=( K2 Poq2 f + ^StY (28) 



12 y/6a 

The FRW evolution laws are regained for (i — |, and t « which corre- 

sponds to high energy densities. Thus we get, for radiation R cx i 1 / 2 , and for dust 
R cx t 2 /3. 

Using the phenomenological limit, a > 1 (MeU) 4 , 29 the inequality, i << 



gives i << 10 30 s. This means that our model is valid for any time in the evolution 
of the universe. The value of q depends on /3. Assuming that the rate of matter 
creation is insignificant during later stages of the evolution of the universe, the 
evolution of scale factor reduces to the FRW relation. 



6. Scalar Field Dynamics 

The effect of matter creation may be translated into the dynamics of a minimally 
coupled scalar field. The effect of negative pressure due to the matter creation is 
equivalent to a suitable scalar field potential. Simple models of a super accelerated 
universe include a phantom field, i.e., a minimally coupled scalar field with negative 
kinetic energy. 30-32 It has also been proposed that early inflation and late time 
acceleration of the universe can be unified in a single theory based on a phantom 
field. 33 Phantom type of matter may also arise from a bulk viscous stress due to 
particle production. 34 

We consider particle production from a phantom field, and derive the corre- 
sponding potential. The field equation including the phantom field, C is, 

R^—Rg^ = 8ttG (V, + P m ) U^U V + Pg^ - d^Cd v C + ^g^g af} d a Cd p C^j (29) 

The RHS can be written as, 8ttG (r^ md + T^ antom ^j , then the sum = 

T^ uld + XP^ antom is conserved by the Bianchi identity but the individual terms arc 
not. This was interpreted in the Steady State theory 35 as creation of matter by the 
creation field in order to maintain a constant matter density. 
The conservation equation, assuming Friedmann metric, is, 

p + 3H(p + p) = (30) 
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We assume the total energy and pressure consist of matter and phantom field 
components as, 

P = Pm+P<f> (31) 

Pm + P<j> + 3H(p m + + p m + p$) = (32) 

Pm + (^p m + p m + -^f>4> + P<t> +P^j = ( 33 ) 

To obtain a phenomenological classical description of the essentially quantum 
creation process, we assume that the creation pressure is proportional to the energy 
density and pressure of the scalar field as, 

Pc = 7^P4> + P4>+P4, (34) 

Now the Eq. (33) becomes, 

p m + 3H (p m + p m + p c ) = (35) 
The expression for creation pressure, from Eq. (11) is, 

Pc = -(Pm+Pm)-^jj = -{Pm+Pm)P (36) 

For perfect fluid matter, p m = (7 — l)p m 



The solution is, 



R 

Pm + 3-7Pm(l-/3) = (37) 



PmO / ocA 
Pm = (38) 



where q = 37 (1 — /3) For phantom fields, 36 ' 30 the energy density, 



P* = -\$ i + V(<t>) (39) 



and pressure, 



P* = -\P-V{<I>), (40) 
where V is the potential of phantom field. The scalar field equations are, 



|(i? 3 0) = O (42) 



The solution is, 



<Po 
i? 3 



(43) 
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Using this solution and Eqs. (24), (25) and (26) in Eq. (23) we get the potential as, 

/ • \ 9/3 

q Ri q Uo J ' 



The modified Friedmann equation, Eq. (1), for very high energy densities becomes, 

R 2 ~ 6a 1 



R 2 K 2 

H 2 = % = 7-J (45) 



Using Eq. (31) for p, 

R ±K 



R \/6a 
We can write this equation as, 



(Pm + P4>) (46) 



R A B 

+ 7^ = (47) 



(48) 



R Ri R 6 

where 

Vfo V Q ) qVfo ' 2\/6^ 

Defining, R q — X where q — — (3), the above equation can be integrated as, 

This integral can be solved for a given value of q. With (3 = 1/2, (the value 
which gives the FRW law for the RSII case), and 7 = 4/3 (radiation), q becomes 2, 
thus the solution of Eq. (49) is, 

X Bq arctan/i ( J ^X ) 

'-"'Ai iAqf r>% (50) 

Neglecting the second term in comparison with the first, and t = we get, 

R = 2A t 1 ' 2 = ^£L t 1 ' 2 (51) 



Substituting in Eq. (43) and solving for we get, 

t= m^L 6 ^ (52) 

(3 K PmOl) <P 

Using Eqs. (44), (51) and (52), we get the expression for the potential in terms of 
phantom field as, 

9/?K 4 (p ro07 ) 5 , 2 . . 

V = 32^1 (53) 

This potential of the phantom field gives the same evolution for the sale factor 
as produced by the adiabatic matter creation, in the high density regime. In this 
sense the phantom field model is equivalent to the adiabatic matter creation model 
discussed previously in section 4. 



4, 2008 7:30 WSPC/INSTRUCTION FILE matter 



9 



7. Gauss-Bonnet brane model 

Gauss-Bonnet (GB) combination arises as the leading order for quantum corrections 
in the heterotic string effective action. In five dimensions it represents the unique 
combination of curvature invariants that leads to second-order field equations, linear 
in the second derivatives in the metric tensor, and it is ghost-free. 37-44 The GB 
correction to the Einstcin-Hilbcrt action will disrupt mildly the degeneracy between 
tensor and scalar perturbations. 45 Recently there is a renewed interest upon this 
scenario. Here we consider the effect of adiabatic matter creation in this scenario. 

In the RSII scenario with a GB term, the effective Fricdmann equation describ- 
ing the motion of brane containing general perfect fluid may be derived from a 
generalization of the Birkhoff's theorem. 46 The 5-dimensional bulk action for the 
GB braneworld scenario is given by, 47 

S = 7T2 / ^x^f—g (r - 2A< 5 > + a [GB]) + S dM + S m (54) 

where [GB] = (R 2 - 4:R AB R AB + R ABCD R ABCD ) is the Gauss-Bonnet term, a > 
with dimension (mass) 2 is the Gauss-Bonnet coupling, A' 5 ) < is the bulk 
cosmological constant and n\ — SttM^ 3 determines the 5d Planck scale. The Sqm 
is the boundary term that is required to cancel normal derivatives of the metric 
tensor which arises when varying the action with respect to the metric, 48 and S m 
is the action for matter confined on the brane. We assume a Z2 symmetry across 
the brane. We will consider the case that a perfect fluid matter source with density 
p is confined to a brane with tension a. 

A new constant is defined, b = (l + |aA( 5 )) . Since A< 5 ) < 0, there is an upper 
bound on a, a < - 1 ^ (5) | = a u . In order that the standard Fricdmann equation be 

2 

recovered at sufficiently low energy scales (p « a), k 2 = = > where M4 is 
the 4d Planck scale. Introducing a dimensionless variable \i the modified Fricdmann 
equations are 



ab 1 / 3 
3 an 2 



1/2 



sinhx (55) 



and, 



H 2 = ±- 

4a 



6 V3 cosh f % ) _ ! 



(56) 



where a is the brane tension. In the high energy approximation, p >> a, we can 
write it as, 

H 2 = Ap 2 ' 3 (57) 

In order for that the 4d effective cosmological constant vanishes, the brane tension 
should satisfy a — § 1 ~^, 2 . We can see that a ss |A^ 5 ) \/n 2 f° r anv < a < a u . 



4, 2008 7:30 WSPC/INSTRUCTION FILE matter 



10 



7.1. Mater Creation on GB Model 




Following a similar calculation as in RS II model in section 5, with a substitution 
X = Rh in Eq. (57), we get 

X* = ttZfAx (58) 
The solution with the boundary condition, X = R q = when t = 0, gives 

(59) 

where g = 3 7 (1-3^) and f3 = ^. 

This evolution is for high energy densities t << 1 , since we considered the 
modified Fricdmann equation for that regime. The case of (3 = 1/2 , which was 
considered in RSII model, gives for radiation and dust as R oc i 3 / 2 and R <x t 2 
respectively. These are modified evolution equations compared to standard cases 
without matter creation. 



8. DGP Brane Model 



An alternative scenario to the RS II brane model was proposed by Dvali, Gabadadzc 
and Porrati (DGP) . 49-51 In this model the extra dimension is infinitely large. The 
key feature is the presence of a four-dimensional curvature scalar in the action. The 
DGP model predicts that 4D Newtonian gravity on a brane world is regained at 
distances shorter than a given crossover scale r c (high energy limit), whereas 5D 
effects become manifest above that scale (low energy limit). The interesting feature 
is that the effective cosmological constant on the brane can be extremely reduced 
in contrast to the case of the Randall-Sundrum model even if a bulk cosmological 
constant and a brane tension are not fine-tuned. Also, the model can explain late- 
time acceleration without having to invoke a cosmological constant or quintessential 
matter. 52 Here we consider a generalized DGP brane model of Sahni and Shtanov. ? 
The action of the generalized DGP model is/ 



where, 



" ,t = / 1 ^[^ ,5,U ™ 

d 5 



, and, 



Sbr 



M 



— K± + Lbrane {9a/3,X) 



(61) 

(62) 
(63) 
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Here k\ is the 5-dimensional gravitational constant, and Lm are the 5- 
dimensional curvature scalar and the matter Lagrangian in the bulk, respectively. 

(fi = 0,1,2,3) are the induced 4-dimensional coordinates on the brane, K ± is 
the trace of extrinsic curvature on cither side of the brane and Lb rane (ffa/3, x) is 
the effective 4-dimensional Lagrangian, which is given by a generic functional of the 
brane metric g a p and matter fields \ 011 the brane. 

The brane Lagrangian consists of the following terms, 



Lbrane = 7T"2 ^ ~~ ° + L m (64) 

where a is the brane tension. 

The 5-D bulk space includes only a cosmological constant A' 5 ) . It is a generalized 
version of Dvali-Gabadadze-Porrati model, which is obtained by setting a = as 
well as A' 5 ) = 0. From the field equation induced on the brane, the effective 4- 
dimensional cosmological constant on the brane is, A = \ (A' 5 ' + 

Considering a Friedmann-Robertson- Walker metric on the brane, the Friedmann 
equation is found to be, 52 

H 2 + ^ = ^[p + p (l + 6A(p,R))] (65) 
where S is either +1 or -1, and A is given by, 



A = 



2 ( 2r) ( n e 



Ai + [p 



(66) 



Po V 4R\ 

with A = y/l - 2?7^, 7] = (0 < n < 1) , and p = m 4 . + 6mf k\. There 

arc two mass scales, m a = er 1 / 4 and m 5 = k^ 2 ^ 3 , in this model. 

At high energy limit, p >> po, the modified Friedmann equation Eq. (65) reduces 

to 



2?7 / e 



For simplicity, the curvature term and dark radiation term can be omitted. Then 
the Eq. (67) can be further simplified to 

# 2 = f [p + S^hfpVp)- (68) 
We will use this evolution law to investigate the effects of adiabatic matter creation. 

8.1. Matter Creation on DGP Model 

We assume the expression , ip = /3 3n H, for matter creation. Defining, q = 37 (1— /?), 
and inserting the solution, p = po R~ q , of the conservation equation (with modified 
creation pressure term) in Eq. (68), we get 

# _ Po4 R - q _ s 4V^£o R - q/2 = (69) 

R 2 3 3 y ' 
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Defining R q / 2 = X, we get 

it x 



R X q/2 
Substituting this in the above equation, 



(70) 



o r =° ( 71 ) 



XV/4 3 X 2 3 X 

This equation can be written as, 

X 2 = AX + B (72) 
Where A = and B = ^p-^f.. The solution of Eq. (72) is 



t-t = ^X + l (73) 

With the boundary condition, X = R q = 0, when t = 0, the constant £ becomes 
_27b_ Therefore we get, 



- t 

3nH- 



where q = 37 (l - ^) , 5 = ±1 and /3 = 

1 /r 

For P = ^ (special solution in the case of RSII), and for ( « ™ , which 
corresponds to high energy densities, we get, 

for radiation, Rent (75) 

for dust, i?cxi 4/3 . (76) 

These are new evolution equations different from FRW laws. For t » (later 
stage) the rate of matter creation is assumed to be insignificant, thus the evolution 
of scale factor is unaffected. 



9. Discussion 

The adiabatic matter creation rate is found to affect the evolution of scale factor 
and energy density of the universe, in all the cases we have considered. The modified 
Fricdmann equation on brane world shows that the universe evolves slowly in the 
very early times. Thus there is a transition of the universe evolution between an 
early high-energy regime and a late low energy regime. The negative pressure due 
to the matter creation process can explain this transition. For the RSII model the 
standard FRW evolution law can be obtained for a matter creation rate (3 = 1/2. 
This value is consistent with the observed data. We have derived a scalar field 
equivalent formulation using a phantom field. A potential is derived which would 
give the dynamics same as that produced by adiabatic matter creation. In the case 
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of Gauss-Bonnet and warped DGP models there are modifications in evolution of 
the universe with matter creation. For later stages of the evolution of the universe, 
the rate of matter creation is insignificant, thus the evolution of scale factor is 
unchanged. 
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